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APPLICATION  OF  THE  DUTCH  METHOD  FOR  ESTIMATING 
STORM-INDUCED  DUNE  EROSION 

PART  I:  INTRODUCTION 

1.  When  a  storm  generates  large  wave  and  surge  conditions,  changes  in 
the  beach  topography  and  nearshore  bathymetry  occur,  often  with  significant 
beach  and  dune  erosion.  Until  recently,  little  quantitative  information  on 
the  prediction  of  beach  and  dune  erosion  existed  to  aid  coastal  engineers. 

The  US  Army  Engineer  Coastal  Engineering  Research  Center's  (CERC)  Shore 
Protection  Manual  (1977)  gives  beach  erosion  quantities  based  only  upon  a 
qualitative  description  of  a  storm. 

2.  Based  on  the  earlier  work  of  Edelman  (1968)  and  comprehensive  model 
testing,  Vellinga  (1983)  reported  on  an  empirical  method  developed  at  the 
Delft  Hydraulics  Laboratory,  Emmeloord,  Netherlands,  for  estimating  dune  ero¬ 
sion  for  a  given  design  storm.  Though  the  method  was  designed  for  the  Dutch 
coast,  it  provides  reasonable  results  for  the  available  data  from  events  along 
the  Gulf  and  Atlantic  coasts  of  the  United  States,  making  It  a  useful  tool  for 
the  design  of  dunes  and  beach  fill  projects  and  for  identifying  beaches  vul¬ 
nerable  to  storm  damage.  The  purpose  of  this  report  is  to  describe  the  method 
and  identify  Its  limitations.  The  report  includes  an  example  problem  using 
the  method  and  an  interactive  computer  program  (FORTRAN  77)  for  applying  it. 

3.  The  three  basic  assumptions  underlying  the  method  are  that  (a)  the 
shape  of  the  poststorm  profile  is  in  equilibrium  (the  procedure  does  not  in¬ 
clude  time  dependency),  (b)  transport  is  in  the  offshore  direction,  and  (c) 
the  amount  of  material  eroded  must  equal  the  amount  deposited.  There  is  no 
provision  for  handling  longshore  gains  or  losses  to  the  profile.  The  con¬ 
cept  is  illustrated  in  Figure  1. 
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Figure  1.  Definition  sketch  of  prestonn  and  predicted  poststorm 
profiles  showing  the  resulting  erosion  and  accretion  zones 

(Vellinga  1983) 


PART  II:  METHOD  OF  ANALYSIS 


4.  The  method  requires  information  about  the  shape  of  the  prestorm  pro¬ 
file  (x^,  y^,  see  Figure  1);  the  sediment  fall  velocity,  w  ;  the  significant 

deepwater  wave  height,  H  ;  and  the  peak  storm  surge  level,  S  .*  The  pre- 

o  s 

storm  profile  should  consist  of  field  data  taken  offshore  to  a  depth  at  least 

equal  to  0.75  H  below  S  (paragraph  6).  If  actual  offshore  data  are  un- 
o  s 

available,  it  may  be  possible  to  estimate  the  offshore  portion  of  the  profile 
using  hydrographic  charts  or  a  predicted  equilibrium  shape  (Dean  1977,  Everts 
1978).  This  will,  however,  decrease  the  accuracy  of  the  method.  The  fall 
velocity,  w  ,  should  correspond  to  the  median  sediment  diameter,  D^^  ,  which 

is  representative  of  the  section  of  the  profile  that  is  expected  to  erode. 
Typically,  D,.^  should  be  computed  based  on  a  composite  of  samples  from  the 
beach  and  dune  zones.  The  fall  velocity  can  be  graphically  determined  from 
Figure  2  for  a  given  water  temperature  and  D^^  . 

5.  The  method  is  applicable  within  the  varying  range  of  wave  steepness, 
0.02  <  H0S/L0  1  0.04  ,  where  Lq  is  the  significant  deepwater  wave  length, 
and  assumes  a  storm  surge  hydrograph  similar  to  those  occurring  in  the  North 
Sea  (Figure  3),  with  S  defined  as  the  summation  of  the  astronomical  high 
tide  and  the  storm-induced  water  level  relative  to  mean  sea  level  (msl) .  This 
hydrograph  is  characterized  by  its  height  and  short  duration,  and  is  similar 
to  that  produced  by  tropical  storms  along  the  coast  of  the  United  States.  A 
procedure  recommended  by  Vellinga  (1983)  for  modifying  the  results  for  longer 
duration  storm  surges  is  discussed  in  paragraph  14.  Ideally,  both  S  and  the 
wave  height  should  be  measured  outside  the  breaker  zone  and  the  corresponding 
deepwater  wave  height,  HQg  ,  should  be  computed.  Estimates  of  Hqs  and  S 
could  also  be  obtained  from  historical  data,  statistical  analysis  (Van 

de  Graaff  1983),  or  numerical  modeling  of  the  storm  conditions. 


*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  defined 
in  the  Notation  (Appendix  B) . 
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Figure  2.  Fall  velocity,  w  ,  as  a  function  of  the  median  particle 
size,  Djq  ,  and  the  water  temperature,  t  .  Curves  obtained  from 
selected  ratios  of  saltwater  to  freshwater  density,  saltwater  to 
freshwater  viscosity,  and  a  sediment  density  characteristic  of 
quartz  sand  (Hallermeier  1981) 
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Figure  3.  Typical  time  history  of  North  Sea  storm  surge,  with 
a  return  period  of  10,000  years  (Vellinga  1983) 

6.  Once  H  ,  S  ,  and  w  have  been  determined,  the  shape  of  the 
os 

poststorm  profile  can  be  obtained  using  the  following  equation: 


where  x  ,  y  ,  an$  H  are  in  meters  and  w  is  in  m/sec.  The  profile 

o  s 

defined  by  Equation  1  terminates  offshore  at 


which,  substituting  into  Equation  1,  yields 


os 


y  =  0.75  H 


(3) 


Constant  slopes  are  used  to  terminate  the  profile  at  the  shoreward  end  (m  at 
x  =  y  =  0)  and  seaward  end  (m2)  of  Equation  1  and  are  defined  as 

m:  =  -1:1  (4) 


m2  =  -1:12.5  (5) 

Use  of  this  steep  poststorm  dune  face  is  consistent  with  field  observa¬ 

tions,  and  m2  is  arbitrarily  fixed  in  agreement  with  model  tests. 

7.  Once  the  poststorm  profile  is  obtained,  its  position  relative  to  the 
prestorm  profile  is  determined  by  horizontally  shifting  the  origin  (x  =  0  , 
y  =  0)  along  S  until  the  erosion  area  equals  the  deposition  area.  Then 

A  ,  the  gross  area  change  A  (equal  to  2A  since  erosion  is  equal  to 
e  g  e 

deposition),  the  erosion  area  above  the  surge  level  Ag  ,  and  the  horizontal 
dune  recession  R  can  be  found  using  either  the  graphical  or  computational 
means  described  in  Part  VI. 
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PART  VIII:  SUMMARY 


22.  An  empirical  model,  developed  using  data  from  the  field  and  from  an 
extensive  series  of  two-dimensional  large-  and  small-scale  movable  bed  model 
tests,  has  been  proposed  by  the  Dutch  to  estimate  dune  erosion  for  large  storm 
surges.  The  model  has  been  evaluated  using  less  extreme  field  data  from  the 
United  States  Atlantic  and  Gulf  coasts  and  should  provide  useful  estimates  of 
storm  erosion.  However,  the  user  should  be  aware  of  the  method's  inherent 
limitations  (Part  V)  and  assumptions  (Part  I)  when  applying  the  model  to 
site-specific  conditions. 
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a  result  Identical  to  the  one  for  the  hurricane.  Consequently,  a  dune  crest 
width  greater  than  27.5  m  is  adequate  for  either  the  hurricane  or  storm 
conditions . 


b.  Given  for  northeaster: 


S  «  3.0  m 

H  =  8.0  m 
os 

D50  •»  0.35  mm 
t  =  2°C 

duration  *=  9  hours 

Solution:  Since  the  width  of  the  dune  crest  is  unknown,  extend  it  in  the 

shoreward  direction  at  a  constant  elevation  of  6  m.  From  the  computer  program 

3 

(or  graphically) ,  we  find  for  the  hurricane  that  A  =  -155.5  m  /m,  A  *  the 

e  s 

erosion  quantity  above  4  m  as  -79.4,  and  R  M  -85.5  m  ,  where  R  is  defined 
ar  the  poststorm  position  of  the  profile/surge  level  intercept.  Then  assuming 
the  1:1  slope  of  the  dune  face,  the  position  of  the  dune  crest  (y  «  +6  m)  can 
be  determined  as 


-85.5  -  (6  -  4)  -  -87.5  m 

Since  the  prestorm  dune  crest  position  is  at  -50  m,  a  crest  width  of  27.5  m  is 
required.  Note  that  the  value  of  Ag  for  the  hurricane  did  not  have  be  dura¬ 
tion  adjusted  since  the  storm  duration  (time  the  surge  is  within  1  m  of  peak 
surge,  S  )  is  equal  to  the  North  Sea  storm  duration  of  5  hours. 

21.  For  the  northeaster,  the  following  erosion  estimates  are 

3  3 

A^  =  -145.6  m  /m  ,  As  *  the  erosion  quantity  above  3  m  *  87.2  m  /m,  and 

R  »  -72.8  m.  In  following  a  treatment  similar  to  that  for  the  hurricane,  we 

find  the  dune  crest  position  to  be 


-72.8  -  (6  -  3)  -  -75.8  m 

In  order  to  adjust  the  volume  and  R  to  account  for  the  9-hour  storm  duration 
of  the  northeaster,  use  the  procedure  outlined  in  paragraph  14.  Assuming 
(conservatively)  a  10  percent  additional  change  in  volume  per  additional  hour, 
then 


10(9  -  5)  a  40  percent 
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Figure  8.  Equilibrium  profile  at  Duck,  N.  C. ,  with  the  predicted 
profile  computed  from  the  storm  conditions  listed 


Example  2 

20.  Using  the  prestorm  profile  of  example  1  and  the  computer  program  in 
Appendix  A,  find  the  maximum  Ag  and  the  width  of  the  dune  crest  required  to 
prevent  dune  failure  for  the  following  conditions  (this  could  also  be  solved 
graphically) : 

a.  Given  for  hurricane: 

S  =  4 . 0  m 

H  *  9.0  m 
os 

D50  =  0.35  mm 
t  =  27°C 


duration  *  5  hours 


PART  VII:  EXAMPLE  PROBLEMS 


Example  1 

19.  Figure  8  shows  Che  predicted  erosion  caused  by  a  simulated  storm 
passing  CERC's  Field  Research  Facility  in  Duck,  N.  C.  The  storm  produced 
7-m  deepwater  waves  and  a  2.5-m  surge.  Using  the  graphical  technique,  compute 
the  erosion  which  will  occur  if  the  surge  had  been  2.7m.  (The  computational 
solution  to  the  problem  Shown  in  Figure  8  is  included  in  Appendix  A.) 


Given:  H 


D 


os 

S 

50 
t 


7.0  m 
2.7  m 
0. 35  mm 
25°C 


Prestorm  Profile  (m) 


i 

i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Xi 

-500 

-no 

-80 

-60 

-50 

-40 

0 

40 

200 

1200 

3 

3.5 

4 

6 

6 

3 

0 

-2 

-5 

-10 

Solution: 


a. 


d. 


e. 


Given 
w 


D, 


50 

0.0474  m/s 


0.35  mm 


25°C  ,  use  Figure  2  to  determine 


Plot  the  prestorm  profile  on  graph  paper.  Compute  the 
poststorm  profile  using  Equations  1-5  and  plot  on  tracing 
paper.  (Or,  since  w  and  Hqs  are  the  same  as  in  Figure  8, 

both  profiles  can  be  traced,  as  poststorm  shape  is  independent 
of  S  .) 

Position  poststorm  profile  at  2.7  m  on  the  prestorm  profile  and 
slide  it  horizontally  until  the  erosion  and  accretion  areas 
between  the  two  profiles  balance. 

3 

Once  the  areas  balance,  compute  A  =  -78  m  /m  and 
3  6 

Ag  *  -44  m  /m  by  planimeter  or  by  counting  graph  paper 
squares . 

Once  plotted,  the  shift  of  the  surge  level  contour  can  also  be 
determined  as  -19.9  m. 
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PART  VI:  GRAPHICAL  AND  COMPUTATIONAL  PROCEDURES 


17.  The  simplicity  of  the  model  permits  a  graphical  solution  of  the 
method.  An  example  is  given  in  Part  VII.  The  procedure  consists  of 

a.  For  given  values  of  H  and  w  (obtained  from  Figure  2), 

construct  a  template  of  the  predicted  poststorm  profile  using 
Equations  1-5. 

b.  Position  the  predicted  dune  toe  (x  =  y  =  0)  of  the  template 
at  the  surge  level  intercept  of  the  prestorm  profile.  Slide 
the  template  horizontally  shoreward  until  the  erosion  area 
appears  to  equal  the  accretion  area. 

c.  Using  a  planimeter  (or  other  method),  find  the  erosion  (-)  and 
accretion  (+)  areas.  Repeat  the  procedure,  sliding  the  tem¬ 
plate  right  or  left  until  the  areas  balance  within  acceptable 
limits  (e.g.  erosion  area  =  accretion  area  ±  5  percent). 

18.  An  interactive  computer  program  written  in  FORTRAN  77  is  included 
in  Appendix  A  along  with  instructions  and  a  sample  run.  The  computer  program 
has  several  features  to  aid  the  interactive  user  including 

a.  Computation  of  w  ,  when  the  user  supplies  D  and  t 
(Figure  2). 

b.  Use  of  English  or  metric  system  of  measurement. 

c.  Multiple  interactive  runs  on  the  same  profile,  changing  S  , 

H  ,  and  w  as  desired, 
os 

d.  Output  of  erosion  quantities  to  the  user's  terminal  for 
immediate  use. 

e.  Error  detection  of  input  data,  failure  of  convergence  on  a 
solution,  or  computation  of  nonvalid  solution. 


i 


model  solution  will  be  valid  if  the  elevation  behind  the  dune  is  signifi¬ 
cantly  higher  (1  to  2m)  than  the  surge  elevation,  causing  little  or  no  over¬ 
wash.  Determination  of  the  dune  failure  point  for  existing  dune  cross  sec¬ 
tions  and  design  dune  cross  sections  under  various  storm  intensities  can  also 
be  useful  in  designing  protective  dunes. 

16.  Another  possible  limitation  which  must  be  considered  is  the  shape 
of  the  natural  and  predicted  profiles.  The  method  should  work  best  for  pro¬ 
files  that  are  similar  to  the  Dutch  ones  used  to  develop  the  model  (Figure  1). 
It  may  not  work  on  profiles  which  deviate  significantly  from  the  typical  Dutch 
profile.  This  problem  was  encountered  first  on  a  beach  with  a  very  flat  off¬ 
shore  zone  which  covered  a  layer  of  peat.  For  this  particular  locality,  the 
profile  shape  predicted  by  Equation  1  was  too  steep  and  always  fell  below  the 
actual  profile  in  the  offshore  region  where  deposition  (not  erosion)  was  ex¬ 
pected.  Similar  problems  may  occur  when  the  procedure  is  applied  to  a  design 
beach,  berm,  or  dune  cross  section  which  may  be  significantly  different  from 
the  natural  profile  shape. 


PART  V:  LIMITATIONS  OF  THE  METHOD 


13.  The  method  assumes  that  all  transport  is  in  an  offshore  direction, 
that  there  is  sufficient  material  to  erode,  and  that  the  poststorm  shape  is  in 
equilibrium.  While  these  assumptions  are  probably  acceptable  for  a  single 
large  event,  over  the  long  term  (or  where  there  are  gradients  in  longshore 
transport)  they  may  not  be  valid.  For  instance,  as  a  result  of  the  equili¬ 
brium  assumption,  for  any  additional  erosion  to  occur  following  a  storm,  a 
larger  storm  must  occur.  In  reality,  a  minor  storm  or  several  minor  storms 
may  cause  the  same  amount  of  change  as  predicted  for  the  larger  event.  There¬ 
fore,  the  engineer  should  ve  aware  of  long-term  erosion  trends  and  use  the 
method  with  care  in  areas  with  known  gradients  in  long-shore  transport  (near 
inlets,  rivers,  man-made  or  natural  barriers,  etc.).  Similarly,  the  method 
should  not  be  applied  to  beaches  without  dunes  or  beaches  with  seawalls. 

14.  Another  assumption  of  the  model  is  that  the  storm  surge  hydrograph 

should  be  similar  to  that  shown  in  Figure  3.  Although  there  is  reasonable 

agreement  between  the  model  data  and  the  United  States  data  (Figure  5),  the 

model  appears  to  work  best  for  large  surge  events  (i.e.,  large  A  ,  Figure  4) 

s 

since  it  was  derived  from  experimental  and  field  data  with  S  greater  than 
3  m.  The  assumed  duration  of  the  surge  peak  (Figure  3)  is  similar  to  that  of 
hurricanes  occurring  along  the  United  States  coast,  but  can  be  significantly 
shorter  than  that  produced  by  a  "northeaster."  Based  on  large  scale  labora¬ 
tory  tests,  Vellinga  (1983)  suggests  increasing  the  erosion  quantity  Ag  5  to 
10  percent  for  each  additional  hour  of  storm  duration  (maximum  additional 
change  not  to  exceed  50  percent  of  Ag  ) .  Storm  duration  is  defined  as  the 
amount  of  time  the  surge  level  is  within  1  m  of  S  ,  minus  the  duration  of  the 
North  Sea  hydrograph  (Figure  3)  under  the  same  criteria  (5  hours).  Note  that 
this  adjustment  should  be  used  as  an  upper  estimate  for  erosional  losses  and 
may,  in  fact,  be  reduced  by  other  factors  such  as  offshore  formations  (para¬ 
graph  12).  For  example,  although  some  of  the  data  shown  in  Figure  5  should 
have  been  duration  adjusted,  a  better  agreement  was  obtained  using  adjusted 
values . 

15.  Model  results  may  become  questionable  at  the  point  of  "dune  fail¬ 
ure,"  when  the  poststorm  profile  slope  m^  intercepts  the  prestorm  profile 
landward  of  the  dune  crest.  The  possibility  of  overwash  and  shoreward  move¬ 
ment  of  sand  Is  not  taken  into  account  in  the  present  model.  Of  course,  the 

14 


\  V  v 


„■  V  V  vvv'.', 


V 


-120 


£ 

o' 

E 

< 


> 

UJ 


13 

tc 

3 


> 

O 

as 

< 


O 

IT 


-100 


-80 


-60 


-40 


-20 


0.20 


0.25 


0.30  0.35  0.40  0.45 

MEDIAN  SEDIMENT  SI2E.  Dso,  mm 


0.50 


Figure  6 
size  D 


50 


Erosion  above  surge  level  A  as  compared  to  the  median  sediment 
and  the  significant  "deepwatlr"  wave  height 

tions  shown 


H 


os 


for  the  condi- 


Figure  7.  Erosion  above  surge  level  as  compared  to  Ag  surge  level  S  and 
the  significant  deepwater  wave  height  HQg  ,  for  the  conditions  shown 
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PART  IV:  SENSITIVITY  OF  THE  METHOD 


11.  The  explicit  termination  of  the  poststorm  profile  (m^ ,  m^)  and  the 
arbitrary  nature  of  the  prestorm  profile  (dune  shape,  barred  profile,  etc.) 
make  it  difficult  to  quantitatively  demonstrate  how  Ag  (or  Ag  )  changes 
with  changes  in  Hqs  ,  S  ,  or  w  .  Vellinga  (1983)  presents  a  set  of  experi¬ 
ments  showing  the  effect  of  several  factors  on  A  .  He  found  that  A  was 

s  s 

proportional  to  HQg  ,  S  ,  and  the  dune  height,  and  that  the  trends  predicted 

by  the  model  agreed  with  the  experimental  data.  Using  several  profile  shapes 

and  varying  H  ,  S  ,  and  w  ,  qualitative  estimates  on  the  model's  varia- 

os 

bility  were  obtained  using  the  FORTRAN  program  given  in  Appendix  A.  In  general, 
A^  and  Ag  vary  directly  with  HQg  and  S  ,  and  are  inversely  proportional 

to  w  .  These  trends  are  shown  in  Figures  6  and  7  for  a  representative  pro¬ 

file  shape.  Depending  on  the  profile  shape,  small  changes  (10  percent)  in 
HQg  ,  S  ,  or  w  can  result  in  either  large  (greater  than  30  percent)  or 

small  changes  in  A  (and  A  ) .  This  points  to  the  equal  importance  of  H 

e  s  os 

S  ,  and  w  in  the  modeling  of  dune  erosion  (although  HQg  and  S  are  im¬ 
plicitly  coupled).  The  water  temperature  t  ,  required  in  determining  w  , 
can  be  very  important  in  determining  Ag  ,  particularly  from  summer  to  winter 
when  a  change  in  t  from  25°C  to  5°C  can  cause  A^  to  increase  by  70  percent. 

12.  Vellinga  (1983)  also  indicated  that  the  method  is  sensitive  to  the 
presence  of  an  offshore  bar  and  trough.  He  found  from  a  series  of  model  tests 
that  when  offshore  bars  were  present,  the  model  overpredicted  Ag  by  28  per¬ 
cent  compared  with  9  percent  overprediction  for  similar  experiments  without 
offshore  formations.  This  discrepancy  may  result  from  dissipation  of  incident 
wave  energy  by  the  bar.  Because  of  the  model's  sensitivity  to  offshore  fea¬ 
tures,  knowledge  of  the  typical  bar  and  trough  formations  in  a  particular 
study  area  may  help  in  interpreting  the  dune  erosion  estimates. 


of  the  prestorm  and  poststorm  offshore  profile  shape,  wave  height,  storm  surge 
and  sediment  size. 

10.  Since  above  surge  level  changes  (A  )  were  small  and  similar  in  mag- 

s 

nitude  to  the  models’  accuracy  (o  ),  values  of  Amg^  [the  change  above  mean 

s 

sea  level  (msl)  ]  have  been  plotted.  Values  of  Amg^  have  not  been  adjusted 
for  duration  (see  paragraph  14).  Even  for  these  relatively  low  quantities, 
the  model  produces  reasonable  estimates  of  erosion  which  are  within  the  accu¬ 
racy  given  in  Equation  6.  Predicted  poststorm  profile  shapes  were  not  usually 
in  good  agreement  with  actual  poststorm  profile  shapes.  Possible  reasons  for 
this  could  be  the  low  magnitude  of  the  surge  levels,  the  use  of  assumed  off¬ 
shore  data  for  some  points,  and  the  timing  of  the  surveys  (some  of  the  data, 
particularly  data  from  Westhampton,  showed  evidence  of  poststorm  recovery). 

The  data  in  Figure  3  suggest  lower  bounds  for  the  storm  surge  and  wave  height 
of  S  =  1.5  m  and  H  =  3.5  m.  At  values  higher  than  these,  there  is  better 
agreement  between  predicted  and  actual  changes.  At  lower  values,  erosion  is 
not  necessarily  universal,  and  some  profiles  on  a  given  beach  experience 
accretion. 


The  method  also  reasonably  predicted  changes  (S  =  2.4  m)  in  Florida  (Hughes 
and  Chiu  1981)  caused  by  Hurricane  Eloise.* 

9.  To  further  confirm  the  method,  historic  CERC  data  from  United 
States’  beaches  for  three  storms  were  analyzed  and  are  presented  for  a  number 
of  profile  lines  in  Figure  5,  along  with  the  Hurricane  Eloise  data  of  Hughes 
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Figure  5.  Comparison  of  predicted  erosion  A  ^  and 
measured  erosion  for  US  Atlantic  and  Gulf  coast  data 

and  Chiu  (1981).  Profile  data  for  the  two  February  1972  storms  at  Westhamp- 
ton,  N.  Y. ,  were  from  DeWall  (1979)  using  hindcasted  wave  data  and  a  fitted 
offshore  slope.  The  storm  surge  was  measured  by  a  tide  gage  at  Sandy  Hook, 

N.  J.  Data  from  Long  Beach  Island,  N.  J.,  including  surge  data  from  Atlantic 
City,  N.  J.  and  wave  gage  data  from  near  Sandy  Hook,  N.  J.,  were  reported  by 
Birkemeier  (1979).  The  offshore  portion  of  the  profiles  were  fitted  using 
data  from  Everts  (1978)  and  Dean  (1977).  Data  from  Duck,  N.  C.,  were  col¬ 
lected  near  CERC’s  Field  Research  Facility  and  include  detailed  measurements 


*  Personal  communication  from  Pier  Vellinga,  Subject:  "Verification  of  Pre¬ 
dictive  Computational  Model  for  Beach  and  Dune  Erosion  During  Storm  Surges 
Verification  for  Field  Data  of  Dune  Erosion  Caused  by  Hurricane  Eloise  at 
Walton  County  in  Florida,  September  1975,"  February  1983. 


PART  III:  MODEL  EVALUATION 


8.  Because  of  the  general  lack  of  extreme  event  profile  data,  the  model 
is  primarily  based  on  an  extensive  series  of  small-  and  large-scale  model 
tests.  Results  from  these  experiments  compared  with  results  from  the  predic¬ 
tive  model  are  reported  in  Vellinga  (1978,  1982,  1983)  and  are  shown  in  Fig¬ 
ure  4  for  Ag  .  Figure  4  also  includes  field  data  from  a  storm  surge  event 
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Figure  4.  Comparison  of  predicted  dune  erosion  A  and 
measured  dune  erosion  (Vellinga  1983) 

(S  =  3.0  m)  that  occurred  on  the  coast  of  the  Netherlands  in  1976.  The  field 
and  experimental  data  are  in  good  agreement  with  the  predicted  erosion  quanti¬ 
ties,  and  poststorm  surveys  following  a  1953  Dutch  storm  surge  (S  =  3.9  m) 
confirm  the  shape  of  the  predicted  poststorm  profile.  Using  statistical  anal¬ 
ysis  of  the  data,  Vellinga  (1983)  specifies  the  standard  deviation  of  the  pre¬ 
diction  (shown  as  dashed  lines  in  Figure  4)  as 


a  =  (0.10  A  +  20)  m  /m 


(6) 


Explanation  of  Program 


1.  A  listing  of  "DUNE,"  the  FORTRAN  77  dune  erosion  program,  is  in¬ 
cluded  here,  preceded  by  an  interactive  example  run.  This  700-line  program 
has  been  adequately  documented  throughout  with  comment  cards  to  aid  users  in 
modifying  and  understanding  the  program.  The  program  first  requests  interac¬ 
tive  input  of  the  prestorm  profile  shape,  sediment  characteristics,  and  storm 
conditions.  It  also  asks  for  the  entry  of  two  datums:  one  to  define  the 
surge  height  and  one  for  use  in  volume  computations.  Once  all  required  inputs 
have  been  entered,  the  program  determines  the  poststorm  profile  shape,  com¬ 
putes  the  location  of  the  poststorm  profile  using  a  modified  secant  method, 
and  computes  the  area  changes  (Ae,  As).  If  the  duration  of  the  storm  surge 
exceeds  5  hours,  the  program  computes  an  adjusted  estimate  of  the  erosion 
using  the  procedure  described  in  paragraph  14  of  the  main  text.  Program  out¬ 
put  is  returned  to  both  the  user's  terminal  and  a  local  file  (FORTRAN  Unit  14) 
for  later  printing. 

2.  All  interactive  input  is  from  the  user's  terminal  (FORTRAN  Unit  5). 
User  input  in  the  example  follows  the  input  prompt  "?".  This  may  differ  (or 
be  a  blank)  on  different  computer  systems.  Required  input  data  are  entered  in 
free  format  with  multiple  entries  separated  by  blanks  (as  in  the  example)  or 
by  commas.  "Yes"  or  "no"  questions  require  either  a  Y  or  an  N  response. 
Program  execution  may  be  terminated  at  any  time  by  entering  a  carriage  return 
after  the  input  prompt  (signaling  an  End-of-File  on  input). 

3.  The  user  is  first  asked  to  label  his  data  with  the  following: 

a.  Locality  code:  2-character  beach  identifier 

b.  Profile  number:  range  or  line  number 

c.  Survey  number:  repetitive  survey  number 

d.  Date:  the  date  of  the  survey  (YYMMDD) 

The  user  is  then  asked  if  the  computations  and  data  entry  are  to  be  done  in 
English  units  (a  Y  response)  or  metric  (an  N  response) .  Units  must  be  con¬ 
sistent  throughout  with  the  exception  of  the  sediment  size  which  must  be 
entered  in  millimeters. 

4.  The  user  must  next  enter  the  survey  data  which  define  the  prestorm 
profile  shape.  The  points  are  entered,  one  per  line,  with  distance  (X^,)  fol¬ 
lowed  by  elevation  (Y^) .  Up  to  110  points  may  be  entered  and  they  should  be 
entered  in  ascending  order  (the  program  checks  for  this).  Data  should  be 


entered  accurately  since  the  program  has  no  provision  to  interactively  list  or 
correct  entered  data.  Enter  a  carriage  return  to  terminate  survey  data  entry. 

5.  The  following  information  is  requested  to  begin  the  computation  (a 
similar  listing  is  also  available  interactively): 

a.  Surge  datum:  Vertical  datum,  relative  to  the  profile  datum  used 
to  define  the  relative  positioning  of  the  surge  height  in  m  (or 
ft). 

i).  Surge  height:  Height  of  the  storm  surge  above  the  surge  datum 
in  m  (or  ft) . 

c.  Wave  height:  Significant  wave  height  H  in  m  (or  ft). 

—  os 

d.  Change  datum:  Vertical  datum,  relative  to  the  profile  data, 

above  which  the  erosion  quantity,  A  ,  will  be  computed. 

s 

e.  Storm  duration:  The  duration  in  hours  that  the  surge  level  is 
within  1  m  (3.3  ft)  of  the  peak  surge  height.  This  is  optional. 
Enter  a  zero  if  the  duration  is  unknown. 

6.  The  use  of  two  datums  in  the  calculations  allows  considerable  flexi¬ 
bility  in  the  input  data.  The  surge  datum  allows  the  surge  height  to  be  en¬ 
tered  relative  to  a  datum  other  than  the  datum  of  the  profile  data.  For  exam¬ 
ple,  profile  data  are  often  measured  relative  to  mean  sea  level  (msl  *  0) 
while  tide  data  may  be  measured  relative  to  mean  low  water  (mlw  “0).  In  this 
case,  using  for  surge  height  the  actual  tide  data,  enter  for  the  surge  datum 
the  difference  between  msl  and  mlw  (which  should  be  a  negative  value).  The 

change  datum  defines  the  lower  boundary  for  the  computation  of  A  .  Nor- 

s 

mally,  A  is  computed  relative  to  the  surge  level.  (To  do  so,  enter  the 
s 

height  of  the  surge,  relative  to  the  profile  datum.)  It  Is  also  possible  to 
compute  a  modified  value  of  Ag  relative  to  any  elevation.  For  instance,  to 
compute  the  change  above  msl,  enter  a  zero  for  the  change  datum  (if  the  pro¬ 
file  data  are  relative  to  msl) . 

7.  The  sediment  fall  velocity  may  be  entered  directly  or  computed  for  a 
given  sediment  size  and  water  temperature.  The  fall  velocity  should  be  en¬ 
tered  in  m/sec  (or  ft/sec),  the  median  sediment  size  in  mm,  and  the  tempera¬ 
ture  in  degrees  C  (or  degrees  F) . 

8.  Following  data  entry,  the  program  computes  and  displays  at  the 
terminal: 

a.  The  total  erosion  area  (A  )  -  the  volume  of  the  total  eroded 

—  e 

area,  regardless  of  datum. 


A3 


b.  The  above  datum  change  (A  )  -  the  volume  change  above  the 
change  datum. 

c.  The  horizontal  shift  of  the  change  datum  intercept. 

d.  The  position,  relative  to  the  profile  data,  of  the  poststorm 
surge  level  intercept. 

e.  An  adjusted  value  of  Ag  which  accounts  for  storm  duration  (if 
the  duration  exceeded  5  hours) . 

Note  that  in  determining  the  adjusted  value  of  Ag  ,  the  program  does  not 
check  to  see  if  sufficient  additional  material  actually  exists.  Use  the  ad¬ 
justed  amount  only  as  an  indicator  of  the  possible  erosion.  In  addition  to 
the  terminal  output,  a  more  complete  summary  of  the  run  is  written  to  FORTRAN 
Unit  14  for  later  routing  to  a  printer.  The  Unit  14  output  includes  a  listing 
of  all  input  and  output  data  and  a  listing  of  the  coordinates  which  define  the 
predicted  poststorm  profile  shape  (unadjusted  for  the  storm  duration). 

9.  Following  each  run,  the  user  is  asked  if  another  run  is  desired. 
Since  the  profile  data  need  not  be  re-entered,  multiple  runs  can  be  used  to 
determine  the  effect  of  varying  any  of  the  parameters. 

10.  There  are  instances  when  the  program  will  be  unable  to  compute  a 
solution.  These  include  the  following: 

a.  Input  errors  resulting  from  incorrect  entry  of  the  input  data 
or  use  of  illegal  values. 

b.  A  surge  height  that  is  higher  than  the  highest  point  on  the 
profile . 

c.  A  change  datum  that  is  higher  or  lower  than  the  data  defined  by 
the  predicted  profile. 

d.  A  predicted  profile  shape  that  falls  either  above  or  below  the 
prestorm  shape  in  a  zone  where  the  opposite  should  occur.  For 
example,  the  predicted  profile  should  fall  above  the  prestorm 
profile  in  the  offshore  zone  and  below  the  prestorm  on  the 
beach.  If  this  does  not  occur,  the  procedure  will  not  work. 

e.  If  the  prestorm  profile  has  insufficient  beach  or  dune  width 
for  the  procedure,  the  predicted  profile  will  fall  landward  of 
all  the  prestorm  data  points  and  an  incorrect  solution  will 
result . 

The  program  prints  an  error  message  when  each  of  these  (or  variations  of  them) 
are  encountered.  When  possible,  the  program  will  still  compute  results  but 
the  user  should  use  them  with  care.  The  best  way  to  identify  the  effect  of  an 
error  is  to  plot  the  prestorm  and  predicted  profile  data  to  see  if  a  reason¬ 
able  solution  was  obtained. 


11.  Should  the  program  fall  to  converge  on  a  solution.  It  may  be  useful 
to  examine  the  sequence  of  Iterations  which  occurred.  These  are  printed  out 
in  reverse  order  (from  the  most  recent  iteration)  on  Unit  14  after  each  run. 
The  top  line  is  the  iterations  required  to  determine  the  horizontal  shift  of 
the  profile  (relative  to  the  datum  of  the  profile  data)  and  the  second  line  is 
the  difference  between  the  erosion  and  accretion  volumes  (which  should  be 
nearly  equal  to  zero  for  the  most  recent  iteration). 

12.  As  written,  the  three  subroutines  to  program  DUNE  are  compatible 
with  the  Interactive  Survey  Reduction  Program  (1SRP)  developed  by  Birkemeier 
(1984).*  With  minor  modifications  to  the  ISRP  calling  and  plotting  routines 
(which  has  been  done),  it  is  possible  to  interactively  generate  and  plot  the 
actual  and  pre-  and  poststorm  profiles,  together  with  the  program's  predicted 
profile  (Figure  8).  In  addition,  ISRP  has  much  more  powerful  survey  data 
entry,  modification,  and  correction  capability  than  program  DUNE.  Corps  of 
Engineers  users  interested  in  either  DUNE  or  ISRP  should  contact  the  Coastal 
Engineering  Research  Center. 


*  References  are  located  at  the  end  of  the  main  text. 


INTERACTIVE  PROGRAM  EXAMPLE 


DUNE  -  A  PROGRAM  TO  ESTIMATE  STORM  INDUCED  DUNE  EROSION. 
WRITTEN  BY  FRANCIS  SARGENT.  CERC.  1983. 

THIS  PROGRAM  USES  A  PROCEDURE  DEVELOPED  BY  MELL INGA 
OF  THE  DELFT  HYDRAULICS  LABORATORY,  THE  NETHERLANDS. 

NOTE  -  A  CARRIAGE  RETURN  AT  ANY  TIME  HILL 
TERMINATE  PROGRAM  EXECUTION 


ENTER  2  CHARACTER  LOCATION  CODE  OF  PRE-STORM  PROFILE 
?  NB 

ENTER  PROFILE  AND  SURVEY  NUMBERS 
?  100  34 

ENTER  SURVEY  DATE  (YYHMDD) 

?  340425 

LENGTH  SCALE  IN  FEET  (ELSE  METERS)?  -  Y  OR  N 
?  N 

ENTER  SURVEY  DATA  -  DISTANCE  THEN  ELEVATION 
ONE  POINT  PER  LINE.  ASCENDING  BY  DISTANCE 

TERMINATE  DATA  ENTRY  WITH  A  CARRIAGE  RETURN 


?  -500  3 
?  -110  3.5 
?  -80  4 
?  -300  4 

LAST  POINT  NOT  ASCENDING.  PLEASE  RE-ENTER  — -  ERROR  MESSAGE  FoR 

?  -60  6  point  sphered  out  of 

?  -50  6  ORDCR 

?  -40  3 
?  0  0 
?  40  -2 
?  200  -5 
?  1200  -10 

7  _ -  CARRIAGE  RETURN  ENTERED  TeRMlNATMc^  DATA  EiTTRV 

DO  YOU  WANT  TO  SEE  A  DESCRIPTION  OF  THE 
INPUT  AND  OUTPUT  VARIABLES7 

?  n  — • - eyOer  a  'Y"  to  <3<rr  a 

ENTER  -  SURGE  DATUM.  SURGE  HEIGHT.  HAVE  HEIGHT.  t>E9CRj pTToN  OF  VARIABLES 

AND  CHANGE  DATUM  ( IN  M  ) 


?  0  n.5  7,0  2.5 

ENTK  HOlWsTORH  SURGE  EXCEEDED  1.50  M  (OR  ENTER  0) 

7  a 

COMPUTE  PARTICLE  FALL  VELOCITY  GIVEN  DIAMETER?  -  Y  OR  N 
7  Y 

ENTER  -  MCDIAN  SAND  SIZE  (MM)  AND  HATER  TEMP.  (DEG  0 
7  .40.10 

FALL  VELOCITY  =  .0474  M  /SEC 


TOTAL  EROSION  (AC)  -  -67.7  M3/  H 

ABOVE  DATUM  CHANGE  (AS)  =  -36.1  M3/  M 

HORIZONTAL  SHIFT  OF  CHANGE  DATUM  =  18.2  M 

POSITION  Or  SURGE  LEVEL  =  -52.4  H 


ESTIMATED  EROSION  ADJUSTED  FOR  DURATION 
(IF  SUFFICIENT  MATERIAL  EXISTS) 

EROSION  ABOVE  CHANGE  DATUM  (AS)  =  TROM  -41.56  TO  -46.90  H3/H 


A6 


I 


OUTPUT  PRINTED  TO  UNIT  14 


DUTCH  STORM  EROSION  PREDICTION 

LOCATION  NE  PROFILE  NO.  100 
SURVEY  NO.  34  DATE  (YYHHDD)  840425 


»*t  INPUT  *»* 


PRE-STGRN  PROFILE 

POINT  DISTANCE  ELEVATION 

NUMBER  M  H 

1  -500.0  3.00 

2  -110.0  3.50 

3  -30.0  4.00 

4  -L0.0  6.00 

5  -50.0  A. 00 

6  -40.0  3.00 

7  0.0  0.00 

8  40.0  -2.00 

q  ?00.0  5.00 

10  1200.0  -10.00 


SURGE  HEIGHT  -  2.50  H 

SURGE  DURATION  =  8.00  HRS 
SIGNIFICANT  HAVE  HEIGHT  =  7.00  H 

MEDIAN  SEDIMENT  DIAMETER  *  .40  MM 

HATER  TEMPERATURE  =10.0  DEG  C 
COMPUTED  FALL  VELOCITY  *  .0474  M  / SEC 
SURGE  DATUM  =  0.00  H 

CHANGE  DATUM  -  2.50  H 


'--.y 


. .  .> 


>  .  -‘.v 


***  RESULTS  »«» 


TOTAL  EROSION  !AE)  =  -67.7  M3/  M 

ADOVE  DATUM  CHANGE  (AS)  =  -36.1  M3/  M 

HORIZONTAL  SHIFT  OF  CHANGE  DATUM  =  -13.2  M 

POSITION  Of  SURGE  LEVEL  =  -52.4  M 


ESTIMATED  EROSION  ADJUSTED  FOR  DURATION 
(IF  SUFFICIENT  MATERIAL  EXISTS) 

EROSION  AEOVE  CHANGE  DATUM  (AS)  =  FROM  -41.56  TO  -46.98  M3/M 


POST -STORM  PROFILE 

DISTANCE  ELEVATION 
M  M 


-55.7 

6.00 

-52.4 

2.50 

-50.0 

2.32 

-40.0 

1.70 

-30.0 

1.20 

-20.0 

.77 

-10.0 

.39 

0.0 

.04 

10.0 

-.28 

20.0 

-.59 

30.0 

-.87 

40.0 

-1.14 

50.0 

-1.40 

60.0 

-1.64 

70.0 

-1.38 

80.0 

-2.11 

90.0 

-2.33 

100.0 

-2.54 

110.0 

-2.75 

111.0 

-2.77 

120.2 

-3.50 

THIS  lS  X-CooRt>iN*Tl= 

OF  FOSTCTOKM  SUNE.  CRFST 


SURGE  POSITION  AND  NET  AREA  ITERATIONS  IN 
REVERSE  ORDER  (FROM  LAST  TO  FIRST). 

-52.433  -52.438  -53.761  -49.676 

.000  .045  12.603  -26.298 
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PROGRAM  DUNE 


PROGRAM  DUNE( INPUT . OUTPUT , TAPE5= INPUT . TAPC6-GUTPUT . TAPE14) 

C 

C  PROGRAM  'PMAIN'  USED  TO  ENTER  PROriLE  DATA  AND  CALL 
C  SUBROUTINE  'PREDICT'  WHEN  USED  AS  A  STAND-ALONE  PROGRAM. 

C 

INTEGER  IP . IDATEP , IPROP . NINC. IP i 
CHARACTERS  ANSWER. REPLY 
REAL  P 

COMMON  /DATA/X(liO! ,Y(li0).2(110) .NP0INTS,PQRT!11Q),C0NSL(110).IU 
COMMON  /DLK / IDATE . IT IME . IPROF . ILOC . ISUR . OIST . IUN IT 
COMMON  /PRDCT/P (110. 3),IP. IDATEP. IPROP .NINCJPl 
COMMON  /FALVEL/PD IA , CTEMP . FALLUEL . SURGE . UAVtHT . TEMP . DAREA 
C  DATA  M.F/2HH  .2HFT/ 

C 

C  FIRST  WRITE  PROGRAM  INTRO 
C 

WRITE (6.900) 

900  FORMAT!//. IX. "DUNE  -  A  PROGRAM  TO  ESTIMATE  STORM  INDUCED’. 
tiX.’DUNE  EROSION. ’./.IX."  WRITTEN  BY  FRANCIS  SARGENT.’. 

$"  CERC.  1983.’.// 

♦IX. "THIS  PROGRAM  USES  A  PROCEDURE  DEVELOPED  BY  YELL  INGA*./. 

♦IX. "OF  THE  DELFT  HYDRAULICS  LABORATORY.  THE  NETHERLANDS.’,//. 

SIX'. "NOTE  -  A  CARRIAGE  RETURN  AT  ANY  TIME  WILL’./ 
t.IX. ’TERMINATE  PROGRAM  EXECUTION’.//) 

100  WPITE(i.iOlO) 

READ  (5.1110.  ERR-- 100 .  END=220 )  ILOC 
IP-0 
C 

C  SET  MAXIMUM  NUMBER  OF  SURVEY  POINTS  ALLOWtG.  THIS 
C  MUST  BE  EQUAL  OR  LESS  THAN  THE  ARRAY  SIZE  OF  X  AND  Y  ARRAYS 
C 

NP0INTSS10 
110  WRITER.  1020) 

READ(5.*.ERR= 110. END-220) IPROF. ISUR 
115  WRITE!*. 1025) 

READ (5, « . ERR- I 15, END-220) IDATE 
1Z0  WRITE(6. 1050) 

READ (5. 1 J00.ERR-120. END = 220 ) ANSWER 
IF (ANSWER. EG. 'Y'lTHEN 
I15NIT=2HFT 

ELSE  IF ( ANSWER . EQ . 'N ' ) THEN 

IUN  IT=2HM 

ELSE 

GOTO 120 

END  IF 

140  HRITE(6,1Q40! 

DO  170  I-l.NPOINTS 

145  READ(5.».ERRd60.ENDd80)Y(I).ZU) 

IF( I.EQ.DGQ  TO  170 
C 

C  CHECK  TQ  BE  SURE  DISTANCES  ARE  ASCENDING 

r 

IF(Y!I).GE.Y(I  i))GO  TO  170 
WRITE(6,150) 

150  FORMAT! ’LAST  POINT  NOT  ASCENDING.  PLEASE  RE-ENTER") 

GO  TO  145 
160  WRITER. 165) 

365  FORMAT! ’INPUT  ERROR. PLEASE  RE-ENTER") 

GO  TO  145 
170  CONTINUE 
C 

C  CLEAR  END-OF-FILE  FLAG  BEFORE  CONTINUING. 

C 

130  X7=E0F<5) 

NPOINTS- I -1 
C 

190  CALL  PREDICT 
C 

200  WRITE16.1030) 

RE AD ! 5 . i 1 00 , ERR =200 , END*  220 ) ANSWER 
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2iO 


IF ( ANSWER. EQ. 'Y ' ) THEN 
WRITE (6 . J  060) 

READ ( 3 . 1 1 00 , ERR=21 0 , END=220 ) ANSWER 
IF  I  ANSWER . EQ - ' Y ' ) GOTO 1 00 
IF ( ANSWER . EO . ' N ' ) GOTO 1 90 
GQT0210 
ELSE 

I F ( ANSWER . NE N ') G0T0200 
END  IF 

1010  FORHATCENTER  2  CHARACTER  LOCATION  CODE  OF  PRE-STDRM  PROFILE") 

1O20  FORHATCENTER  PROFILE  AND  SURVEY  NUMBERS" ) 

1025  FORHATCENTER  SURVEY  DATE  (YYHHDD)") 

1030  FORMAT ( "DO  YOU  WANT  ANOTHER  RUN?  -  Y  OR  N") 

1040  FORMAT! "ENTER  SURVEY  DATA  -  DISTANCE  THEN  ELEVATION"./. 

*  "ONE  POINT  PER  LINE.  ASCENDING  BY  DISTANCE", 

I  //. "TERMINATE  DATA  ENTRY  WITH  A  CARRIAGE  RETURN",//) 

1050  FORMAT ["LENGTH  SCALE  IN  FEET  (ELSE  METERS)?  -  Y  OR  N") 

1060  FORMAT! "NEW  PROFILE?  -  Y  OR  N") 

1100  FORMAT (At ) 

1110  FORMAT ( A2 ) 

220  STOP 
END 
C 

CtmKflKimilHHIHHmKIlilliimiHItHIHHHHHHHmtimi 

c 

SUBROUTINE  PREDICT 
C 

C  PREDICT  -  CERCFRF  -  BCS:FTN5  -  F.  E.  SARGENT  -  S30609 

C  PROFILE  GIVEN  A  'PRESTORH'  PROFILE  AND  SEVERAL  BULK  PARAMETERS 

C  OF  THE  STORM  AND  BEACH.  THE  PROGRAM  FOLLOWS  A  METHOD  PRESENTED 
C  BY  VELLINGA  (1983).  THE  POSITIVE  Y-DIRECTION  IS  OUTWARD 
C  NORMAL  AND  THE  POSITIVE  2-DIRECTION  IS  UPWARD.  'PREDICT'  CAN  BE 
C  USED  AS  A  STAND-ALONE  PROGRAM  WITH  ITS  CALLING  PROGRAM  'PMA IN '  OR 

C  AS  A  SUBROUTINE  TO  PROGRAM  '  ISRP '  DEVELOPED  BY  BILL  BIRKEMEIER  AT 

C  THE  CERC  FIELD  RESEARCH  FACILITY.  PROGRAM  CODING  IS  IN  F0RTRAN77. 
C 

REAL  DZ1  DZ2.DY.SAREA  CW.CY .CZ . ZACT 1 . ZACT2 , YACT1113 ) , YACT2 . BAREA 
RE  AL  M 1 .  Il2 .  TARE  A 1 .  TARfc  A2 .  DA&EA  ( 1 3 ) ,  T 1 .  SURGfe .  WAVEHT .  FALLVEL .  DATUM 
RE AL  P . Y IT , ZEND . YEND . C ARE A tSLOPE I . YC 1 , YC2 , ZC 1 . ZC2 . Y I . Z I .AREA ( 3 ) . PR 
REAL  DY2.EAR£A.C.Cl.C2.C3,(lATUM2.Rl.K&.PDIA.TEMP,CTEflP^0S(2>.T2 
REAL  SHIN. WHIN, FHIN.PMIN 

INTEGER  INDEX . INDEX1 . INDEX2 . IP . 3DATEP .FLAG.NINC, IP1.I, J.K.L.M.N.ML 
INTEGER  ZS IGN 1 , ZS IGN2 . IY . YS1 . YS2 . YE1 , YE2 . IRODE . ks IGN1 . AS IGN2 . IPROP 
INTEGER  FRNT.QPOUT  U 
CHARACTER 1 2  U2 . U4 . REPLY  * 1 . ANSWER* 1 
C 

C  LABELED  COMMON  CLOCKS  USED  WITH  PROGRAM  'ISRP'. 

C 

COMMON  /DATA/X!110!,Y(11O).Z(110).NPOINTS,PORT(11O),CONSL(11O).IU 
COMMON  /DLK/IDATE.  IT  [ME.  IPROF,  HOC.  ISUR.DIST,  IUNIT 
COMMON  /PRDCT/P ( 1 10. 3) JP . IDATEP . IPROP . NINC . IP1 
COMMON  /FALVEL/PDIA.CTEHP .FALLVEL. SURGE, WAVEHT. TEMP. BAREA 
C 

C  GLOBAL  CONSTANTS  (NOT  CHANGED  ON  SUBSEQUENT  CALLS  TO  'PREDICT') 

C  NINC  <  OR  --  DIMENSION  OF  PROFILE  ARRAYS  IN  LABELED  COMMON  'DATA' 

C  MAY  OTHERWISE  BE  DEFINED  IN  A  'BLOCK  DATA'  STATEMENT. 

C 

DATA  Kl. K2. Ml, M2. T1.T2, ML/. 471405.. 752266.-1., -.08.. 01.. 01.13/ 

DATA  SHIN, WHIN, PHIN.FMIN.QPOUT/3*. 01 . .001 ,0/ 

NINC=U0 
N=N INC  - 1 
C 

C  SET  FLAGS  FOR  'ISRP'  PLOTTING  ROUTINE  'IPLDT'. 

C 

IP-1 

IDATEP- IDATE 
IPROP -IPROF 
U- IUNIT 
C 

C  USE  APPROPRIATE  UNITS  (FEET  OR  METERS!  FOR  ANALYSIS. 

C 


AIO 


If ( 1UHIT.EG.2HM  ) THEN 

C*l. 

Cl=l. 

C2=0. 

C3=2. 

U2='  M' 

U4-'  C' 

ELSE 
C=.304F 
Cl =.5556 
C2=32. 

C3=54. 

U2='YD' 

U4='  F' 

END  IF 

100  CONTINUE 

IFUPD2.EQ.DC0  TO  102 
C 

C  DETERMINE  IF  AN  EXPLANATION  OF  THE  VARIABLES  IS  DESIRED? 

C 

»*» 

WR ITE 16.1011) 

1011  FORMAT!/, IX. "DO  YOU  WANT  TO  SEE  A  DESCRIPTION  OF  THE'./, 

$1X. "INPUT  AND  OUTPUT  VARIABLES?") 
READ!5.1030.END=410,ERR=100)IREPLY 
IF! IREPLY.EQ. 1 HY ) CALL  EXPLAIN 

ipdm 

102  CONTINUE 
C 

C  REQUEST  ENTRY  OF  SURGE  AND  NAVE  HEIGHT  PARAMETERS 
C 

102  rfRITE(6. 1000)  U 

RE  AD!  5 .  * .  ERR=105 .  END=41 0 )  DATUM ,  SURGE .  HAUEHT .  DATUM2 
C 

C  REQUEST  STORM  DURATION  IN  HOURS  FOR  LATER  ADJUSTMENT 
C  OF  THE  EROSION  AMOUNTS  (IF  DURATION  GT  5  HRS) 

C 

S=SURGE-1./C 

106  URITE!6.!07)S,U 

107  FORMAT! IX, "ENfER  HOURS  STORM  SURGE  EXCEEDED", 

JF&.2.1X.A2.1'  (OR  ENTER  0)") 

READ(5.«.ERR=106.END=410)DUK 

C  SUBTRACT  DURATION  OF  NORTH  SEA  STORM 
C 

DUR=DUR-5. 

IHDUR.LT.O.  )DUR=0. 

c 

C  COMPUTE  FALL  VELOCITY  GIVEN  MEDIAN  PARTICLE  DIAMETER  (IN  MM) 
C  AND  HATER  TEMPERATURE.  OTHERWISE  INPUT  FALL  VELOCITY. 

C 

URJTEit.1020) 

READ (5 . 1030. ERR= 100 . END=410)REPLY 
IFfREPLY.EQ. 'Y'JTHEN 
110  HRITE(6.1010)U4 

READ!5.«,ERRi110,END=410)PDIA.TEMP 

C 

C  CONVERT  TEMPERATURE  TO  DEG  C  (IF  REQUIRED) 

C 

CTEMP=C1*(TEMP-C2) 

IF(CTEMP.GE.0..AND.CTEMP.LE.40,.AND.PDIA.GT.PMIN)G0T0120 

HRITE(A.1Z50) 

GOTOliO 

120  CALL  SPEED 

FALLVEL-FALLVEL/C 

HRITE(4.123)FALLVEL,U 

123  FORMAT!"  FALL  VELOCITY  =  \F7.4.1X,A2.7SEC") 

ELSE  IF ! REPLY. EQ.'N') THEN 
PDJA--0. 

130  HR ITET6 . 1015) 


.y.*  i. 


All 


c 

c 

c 

140 


C 

c 

c 

IPO 


160 


170 


C 

c 

c 


c 

c 

c 

c 


180 

170 

c 


RE AD ! 5 . s . ERR = 1 30 . END*  4 1 0 ) F ALL VEL 

ELSE 

GDT0100 

END  IF 

IF ( FALLOEL . GT . FH IN . AND . UAVEHT . GT . HH IN. AND . SURGE . GT . SH IN) GDT0140 
WRITE(6. 1240) 

GOTOIOO 

INITIALIZE  CONSTANTS  AND  FLAGS. 


ASIGN1 =  -1 
ASIGN2--1 

CH-(rALLVEL+C/0.0268)«.56 

CZ=7.6/(HAVEHT*C> 

CY-CZ**1.28»CH 
ZACT1=DATUM+SURGE 
IF ! Z ( NPO INTS ) . LT . ZACT 1 )G0TQ150 
WRITEI6. 1060) 

ASSIGN  1060  TO  FRMT 

FLAG-2 

G0T0340 

INITIALIZE  FIRST  TWO  'GUESSES'  OF  POST-STORM  PROFILE  POSITION. 

DC  160  I-NPOINTS.l.-l 

IF(Z!  I)  .LT.ZACTDG0T0160 

INDEX- I 

GOT0170 

CONTINUE 

URITEI6.1070) 

ASSIGN  1070  TO  FRNT 

FLAG-2 

GQT0340 

SLOPE=(y(I)-Y(M))/(Z(I)  ZII+D) 

YIT= (ZACT1  Z( I+l ) ) tSLOPE+Y( I+t ) 

YEND=250./(CY*C) 

YACT1 ( 1 ) =Y IT-YEND/10. 

YACT1 (2) -Y IT- YEND/8. 

ZEND=K2+NAVEHT 

ZACT2-ZACT1-ZEND 

CAREA- YEND*ZACT  1  +2 . *YEND/ ( CZ«C ) •*•  ( 7Z .  -2. *K  1*( CY*C* YEND+1 B . ) **  1 . 5} / 
4  (3.+C+C+CZ+CY) 


LOOP  INTERAT ION  STARTS  HERE,  USING  A  SECANT  METHOD  TECHNIQUE. 

DC  270  J-i.ML 
FLAG-0 

YACT2= YACT1 ! J ) +YEND 


LOCATE  LOWER  INDEX  OF  THE  PRESTORM  PROFILE  AT  THE  LEFT 
OF  CURVED  PORTION  OF  PREDICTED  PROFILE. 


END 


DO  590  I* INDEX* 1. 1,-1 
IFIYI I).GT.YACT1(J))G0T0190 
SLOm= ( 2  ( I)-Z(I+1))/(Y(I)-Y(I+1)) 
IRSL0PEI.LE.MUG0T019O 

YC1 = (H1+YACT1 ( J)-5LDPEI*Y ( I+l )+Z( I+l ) -ZACT1 )/ (Nl-SLDPEI) 

IFIYCl .LT.Yi I) )G0T0190 

INDEX1--M 

IF( YC1 .LE. YACT1 ( J) )G0T0180 
YCRYACTKJ) 

FLAG--1 

ZCl-SLOPEI+IYCl-YU+lD+ZII+l) 

GOT0200 

CONTINUE 

WRITEI6.1080) 

ASSIGN  !0C0  TO  HINT 

“LAG-2 

G0T034Q 


LOCATE  UPPER  INDEX  OF  THE  PRESTORM  PROFILE  AT  THE  RIGHT  END 


A12 


C  OF  CURVED  PORTION  Of  PREDICTED  PROFILE. 

C 

200  DO  220  1= INDEX. NPOINTS.l 
IFIYt I).LT.YACT2)GOTG220 
SLOPED  (Z£  I  )-Z(  I  1 ))/('/( I)- Y(I-D) 

IF! SLOPE  I . LE . M2 ) G0TQ220 

YC2--  (M2*  YACT2-SL0PE I  *Y  ( I -1 )  >1  ( I  1 )  -ZFCT2 )  / !  M2-SLQPE I ) 

IF(YC2.GT.Y(I))G0T0220 

INDEX2=I-i 

IT ( YC2 . GE . YACT2 ) G0T021 0 

YC2-YACT2 

FLAG= 1 

210  ZC2=SL0PEI* ( YC2- Y (I  •  1 ) ) +Z (I -1 ) 

GOT0230 
220  CONTINUE 

HRITE16. 1090) 

ASSIGN  1090  TO  FRMT 
FLAG =2 
G0T0340 
C 

C  COMPUTE  AREAS  RELATIVE  TO  ZERO  DATUM  -  PREDICTED  AREA. 

C 

230  T ARE A1 = ( ! ZC 1 +ZACT 1 )*( YACT 1 (J ) -YC 1 ) + 1 ZC2+ZACT2 ) * ( YC2  YACT2 ) ) / 

C 

C  PROFILE  AREA. 

n 

TAREAi-TAREAl+CfiREA 

IF( INDEX2 .G2. INOEXI ) THEN 

TAREA2 = i ZC1 1 Z ( INDEX 1 ) 1 » I Y  < INDEX 1 ) -YC 1 ) 

IF( INDEX2.EQ. INDEX1 >G0T0250 
DO  240  I-JNDEX1.INDEX2-1 
240  TAREA2=TAREA2+«Z(I)+Z(  IH))*(Y(  IH)-Yt  I>) 

210  T ARC A2S i T ARE A2 1 ! Z ( INDEX2  H  ZC2 ) * ( YC2-  Y ( INDEX2 ) ) ) /2 . 

TAR'EA2--(ZCl4ZC2)*<YC2-YCl)/2. 

END  IF 

r* 

C  CHECK  IF  AREAS  LESS  THAN  TOLERANCE  (OR  FLAG  IS  EXECUTED). 

C  EXIT  IF  SO.  OTHERWISE  COMPUTE  NEW  GUESS  AND  RECOMPUTE  AREAS. 

C  ITERATION  PROCEDURE  BASED  ON  STANDARD  SECANT  EXTRAPOLATION 
C  TECHNIQUE  WHEN  NET  AREAS  (DAREA)  ARE  OF  THE  SAME  SIGN. 

C  WHEN  NET  AREAS  ARE  OF  OPPOSITE  SIGN  PROGRAM  SELECTS  TWO 

C  BOUNDING  VALUES  TO  INTERPOLATE  NEW  PROFILE  POSITION. 

C 

DAREA ( j ! - TAREA2 -T AREA 1 
IF(ADS(DAREA(  J) )  .LC.T1IGOT028O 
260  IF i DAREA (J).LT.O) THEN 
IF(ASIGN2.EQ.1)THEN 
IT ( DARE A ( J ) . GT . DAREA ( JZ ) ) J2S  J 
ELSE 
ASJGNZ-I 
J2= J 
END  IF 
ELSE 

IFIASIGNl.EQ.DTIiEN 

IF!DAREA(  J).LT.DAREAUI))J1-J 

ELSE 

ASIGN1-1 

J1  =  J 

END  IF 

END  IF 

IF ( ASIRNi . EQ.l . AND .AS IGN2 . CQ. i ) THEN 

JN  -  J 1 

J0-J2 

ELSE 

JN-J 

JOU  1 

END  IF 

IF U. EQ.l  .OR.  J.EQ.HDGQT0270 

SLOP  E  -■  ( YACT  1  UN )  -  YACT  1<  JO) )/( DAREA  ( JN  >  -D  ARE  A  ( JO ) ) 

YACT  1 1 .1+1 )  -  YACT1  i  JN)  DAREA ( JN)  «SLOPE 


270  CONTINUE 
FLAG-3 
J -Hi. 

C  COMPUTE  ABSOLUTE  AREA  CHANGE  AND  CHANGE  ABOVE  AREA  CHANGE 

C  DATUM.  ADDITIONAL  CONVERGENCE  CHECK  BASED  ON  NET  AREA  CHANGE. 

C  DISCRITIZF  PRE/POST  PROFILES  BETWEEN  THEIR  INTERCEPTS. 

C 

280  P!I.1):YC1 

P(1.2i=ZCl 
P(;,3)=ZC1 

IF ■ ftBCO  YACT1 ( J ) - YC1 ) .LT. . 0001 )P ( 1 . 3) =ZACT 1 
PINING.  D-YC2 
PINING. 2!-ZC2 
PININC  2)zZC2 

IF  I ABS ( YACT2-YC2 ) . LT . . 0001 ) P  t  NINC . 3 ) =ZACT2 
C 

C  COMPUTE  DISTANCES  AND  PRE/POST  STORM  ELEVATIONS. 

C 

L= INDEX 1 

DY-YC2-YCI 

DO  300  1=2. N 

P(I.!)=YCltDY»(I-l)/N 

IF  IP  1 1.  D.L  T.YACT1 1  J)  (THEN 

P! I,3)=ZC!+Kl*(Pi I.1/-YC1) 

ELSE  IF ( P ( I . i ) . LT . YACT2 ) THEN 

P(J.3)=ZACT1-(K1JSQRT(CY«C»!P(I.I)-YACT1(J))+I8.)“2.)/(C2*C) 

ELSE 

P(I.3)=ZACT2:M2*(P(I.1)-YACT2) 

END  IF 

290  IF(YiL).GE.P(I.l))GOT0300 
L=L«-i 
GOT0290 

300  P(I.2)=Z(L-i)KPiI.i)-Y(L  i))*(ZIL)-Z(L-!))/(Y(L)-Y!L  D! 

C  SUM  ABSOLUTE.  NET  AREAS. 

C 

ARFA(1)=0. 

OV-DY/N 

DZ2-PI1.2)  P(i.j) 

Z3IGN2-1 

IF!DZ2.LT.0.)ZSIGN2=-1 

C 

C  ABSOLUTE  AREA. 

C 

DO  310  !=Z.NINC 
DZ1--DZZ 
ZSIGNI-ZSIGN2 
DZ2--PII.Z) -PC  1.3) 

Z5IGN2-1 

IF IDZ2.LT.0. )ZSIGN2=- 1 

IF i ZS  5GNZ . EQ . ZS ION 1 ) THEN 

AREA  ■  1 '  -  AREA  ( 1 '  +ABS  ( ( DZ2-<  DZ1XDY) 

ELSE 

AREA!  1 ) -AREA!  1 )  KDZ2*s2^DZl*«2) « ABS( DY/ (DZ2  DZ1  > ) 

END  IF 

310  CONTINUE 
SAREA=AREAI1)/(Z.»C3) 

K=SAREA»iO.+,5 

8AREA-K/10. 

PR -ABS ! DAREA ( J ) /SAREA ) 

C 

C  IF  CHANGE  DATUM  IS  OUTSIDE  POST  STORM 
C  PROFILE.  REQUEST  A  NEW  DATUN2 
C 

31 1  IF ( BATUM2 . GE . ZCZ . AND . DATUM2 . LE . ZC 1 ) GO  TO  319 
31 Z  WfiITE(6.315'< 

313  FORMAT!"  CHANGE  DATUM  IS  OUTSIDE  THE  '  IMITS"./, 

Of  THE  POST  STO01*  PROFILE.'./. 

*"  ENTFR  NEW  CHANGE  DATUM*) 

READ ( 5  <  ERR- 3 1 2 . END- 4 10 ) DATUM2 


GO  TO  311 


319  DO  320  11=2,3 
AREA ( 1 1 ) = 0 . 

POS( I I— 1 )=0. 

7c inu?- i 

IF  (DZ2.LT. 0.  )ZSIGN2=-i 
IR0DE=0 

DO  320  J=2,NINC 
DZ1 =D22 
ZSIGNi=2SIGN2 
DZ2=P<  r.II)-DATUH2 
ZS3GN2=1 

IF(DZ2.LT.0.)ZSIGN2=-1 

IF ( ZS IGN2 . EQ . ZS IGN1 >  THEN 

IF  ( ZS IGN 1 ,  EG .  1 )  ARE  A  ( 1 1 ) = AREA  { 1 1 )*  ( DZ2+DZ1 )  *DY 

ELSE  IFLZSIGNl .EG . 1 )THEN 

AREA  (II)  =AREA(  1 1 )  +DZH*2*DY/(DZ1-  DZ2) 

IF( IRODE. £0,0) THEN 

PQS( I 1-1 >=P(I- 1 , 1 )+DY»DZl/ (DZ1-DZ2) 

IRODE =1 
ELSE 
END  IF 
ELSE 

AREA (II) » AREA (II) +DZ2**2*DV/ ( DZ2-DZ1 ) 

END  IF 

320  CONTINUE 

EAREA= ( AREA ( 3 ) -AREA (2) >/C3 
DP =P03 ( 2  > -POS ( 1 ) 

C 

C  TOLERANCE,  LOOP  CHECK. 

C 

IFIJ.E0.HUG0T0330 
IF (PR.LE.T25GQT0330 
GOT0260 
330  I?=2 

SAREA=-j , *5AREA 
C 

C  COMPUTE  ADDITIONAL  EROSION  BASED  ON  STORM  DURATION 
C  SKIP  COMPUTATION  IF  DUR  =  0 
C 

IF(DUR.E3.0.)G0  TO  335 
AD J3=EAREA*DUR  * . 05+EAREA 
ADJ4=EAREA»DUR* . 1 +EAREA 
IF( AD J3 . GT . EAREA/2. ) ADJ3=EAREA/ 2 . 

IF ( AD  J4 . GT . EAREA/2 . >  AD J4=EAREA/2 . 

C 

335  IF (FLAG.EQ. 1 )HRITE(6, 1150) 

UR[TE(6,1140)SAR£A.U2,U,EAI?EA.UZ,U.DP.U,YACTI(J).U 

IF(DUR.IjT.0.)HRITE(6,iZ70)ADJ3,ADJ4,U2.U 

IF(FLAG.E0.3)HRITE(6,1I30) 

340  IF ( OPOUT . EO . 0 ) G0TQ360 

350  WRITE(6, 1260) 

RE  AD  (5 . 1 030  .ERR*  350 .  ENDM10)  ANSWER 
IF ( ANSWER .EG. 'N ' ) RETURN 
IT ( ANSWER .EG. 'Y')G0T0360 
G0T0350 
C 

C  WRITE  DATA  INPUT/OUTPUT  TO  UNIT  14. 

C 

360  HR ITE (34.1160) 

WRITE! 14, 1040) ILOC. IPROF, ISUR, IDATE 
HR  I TE  ( 1 4 . 1 1 90 )  U ,  U 
DO  370  1=1 .NPOINTS 
370  WRITEU4, 1230)1, Y(I),Z(I) 

IF  (REPL.Y.EG.  'Y')THEN 

WRJTE(14,1110)SUR0E,U,DUR+5.,NA0EHT,UfPDIA,TEMP.U4,FALLVEL.U, 

?DATUH,U,DATUH2,U 

ELSE 

HR ITE ( 1 4 . 1 1 00 ) SURGE , U , DUR  *5 . , WAVEHT , U , F ALL VEL , U , D ATUM . U , 
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$DATUH2.U 
END  IF 

IF  <  FLAG . NC - 2  >  GOT0380 
WRITE! 14.1120) 

WRITE  U4.FRMT) 

RETURN 

300  WRITE (14.1170) 

HRJTE(i4,1140)SAREA,U2,U,EAREA,U2,U,DP,U,YACTl( J) ,U 
IF(DUR.GT.0)WRIT2!14,li70)A0J3.ADJ4,U2,U 
C 

C  CONFUTE  AND  WRITE  PREDICTED  PROFILE  AT  10  FOQT/METER  INTERVALS. 

C 

YSl-YACTK  J)/10. 

YE  1 =YACT2/ 10 . 

IF  (YS1.E3.0.AND.YC1.LT.0. )G0T0390 
IF( YS1 -GE-O) Y31 =YS1 +1 
390  YS2=YS1*10 

YE2= YE1 » 10 
HRITE(14.1130)U,U 
WRITE(14.1210)YC1.ZC1 
WR ITE! 14.1 210 ) YACTl (J5.ZACT1 
DO  400  IY-YS2.YE2.10 
YI-IY 

ZI=ZACT1-(K1*SQRT(CY*C*(YI -YACTK J)  )+18. )  2. >/(CZ*C> 

400  WRITE! 14. 1210) YI.ZI 

WRITE (14, 1210 )YACT2.ZACT2 
WRITE! 14. 1210)YC2,ZC2 
IFiFLAG.EG.l (WRITE (14.1150) 

WRITE! 14, 1200) 

HRITEU4, 1050)  (YACTl  (J<1-N),H-1.J) 
HRirEil4,i050)(DAREA(J+l-H),H*l. J) 

IF (FLAG. EG. 3) WRITE! 14.1 130)T1 .DAREA ( J ) .T2.PR 
WR  ITE ( 14. 1220) 

1000  FORMAT! "ENTER  -  SURGE  DATUM,  SURGE  HEIGHT.  HAVE  HEIGHT,",/. 

4  8X. "AND  CHANGE  DATUM  !IN",1X.A2, *)",//) 

1010  FORMAT ("ENTER  -  MEDIAN  SAND  SIZE  (MM)  AND  HATER  TEMP.  I DEG*. 
$A2,")") 

1015  FORMAT! "ENTER  -  FALL  VELOCITY") 

1020  FORMAT! "COMPUTE  PARTICLE  FALL  VELOCITY  GIVEN  DIAMETER?". 

4  "  -  Y  OR  N") 

1030  FORMAT(Al) 

1040  FORMAT!"  LOCATION  ",A4,"  PROFILE  NO  "I?,/, 

$  "SURVEY  NO.". 15."  DATE  (YYMMDD)  .16.//) 

1050  FORMAT!  1XJ3F10. 3) 

1060  FORMAT!"  £nD  POINT  OF  PROFILE  IS  ABOVE  SURGE  LEVEL") 

1070  FORMAT!"  SURGE  LEVEL  ABOVE  PROFILE") 

1030  FORMAT!"  PREDICTED  CURVE  INTERCEPT  TO  THE  LEFT  OF  PROFILE") 

1090  FORMAT!"  PREDICTED  CURVE  INTERCEPT  TO  THE  RIGHT  OF  PROFILE") 

1100  FORMAT!//."  SURGE  HEIGHT  =  ".F6.2,1X,A2,/. 

4  ’  SURGE  DURATION  -\F6.2,"  HRS"  /, 

4  "  SIGNIFICANT  WAVE  HEIGHT  =  ".ri.Z.lX.ftZ./, 

$  '  MEDIAN  PARTICLE  FALL  VELOCITY  =  \F6.4,lXfA2,7SEC",/. 

*  "  SURGE  DATUM  =  \r6.2.1X,A2,/. 

4  "  CHANGE  DATUM  =  *.F6. 2. lJ<.A5) 

11.10  FORMAT!//."  SURGE  HEIGHT  =  ".F6.2.1XTA2,/, 

4  ■  SURGE  DURATION  =".F6.2.H  HRS",/. 

*  "  SIGNIFICANT  HAVE  HEIGHT  =  ".fo.Z/X.AZ./, 

4  ’  MEDIAN  SEDIMENT  DIAMETER  *  ’.F6.t."  HH*,/, 

*  "  HATER  TEMPERATURE  =  \F6.1,1X."DEG\AZ./. 

4  "  COMPUTED  FALL  VELOCITY  =  \F6.4,1X,A2,7SEC\/, 

%  "  SURGE  DATUM  *  ".F6.2.1X.AZ./. 

4  "  CHANGE  DATUM  =  •  F6.2.1X.A2) 

1120  FORMAT!"  DATA  INPUT  ERROlS  -  CHANGE  SURGE  HEIGHT.",/. 

4  "  PROFILE  LIMITS  OR  CHECK  COORDINATES  OF  DAtA*") 

1130  FORMAT!"  WARNING  -  SOLUTION  DID  NOT  MEET  CONVERGENCE  CRITERIA"./. 
4  ”  AREA  TOLERANCE  *  \F6.3."  NET  AREA  =  ",F6.3./. 

$  "  RELATIVE  TOLERANCE  =  \F6.3  "  RELATIVE  ERROR  =  ’.F6.3) 

1140  FORMAT!///  TOTAL  EROSION  \ AE >  =  n.F8.1.tX.AZ,"3/".lX,A2./, 

J  "  ABOVE  DATUM  CHANGE  (AS)  *  ".F8.1 .1X.A2/3/*. 

41X.A2./. 

4  "  HORIZONTAL  SHIFT  OF  CHANGE  DATUM  •  *,FB.1.1X,AZ./. 


t  "  POSITION  OF  SURGE  LEVEL  =  ".F8.1.1X.A2.//) 

1150  FORMAT ( "  WARNING  -  PORTION  OF  PREDICTED  PROFILE  ABOVE/BELOH*./, 
t  "  PRESTORM  PROFILE  WHERE  REVERSE  TREND  SHOULD  OCCUR*) 

1160  FORMAT  I //,5X, "DUTCH  STORM  EROSION  PREDICTION".//) 

1170  FORMAT! ////.i6X.’m  RESULTS  m*.///) 

1180  FORMAT!"  POST-STORM  PROFILE  *.//. 

$  ■  DISTANCE  ELEVATION”./. 7X.A2.9X.A2./) 

1190  FORMAT 1///.10X.’***  INPUT  ***'.,'//, 
r  PRE-STORM  PROFILE  *.//, 

i  ’POINT  DISTANCE  ELEVAT ION” , /  *  NUMBER* , 7X . A2 . 9X , A2 , / ) 

1200  FORMAT!//. "SURGE  POSITION  AND  NET  AREA  ITERATIONS  IN*,/, 

*  "REVERSE  ORDER  (FROM  LAST  TO  FIRST).*,/) 

1210  F0RMAT(1X,F9, 1.2X.F10.2) 

1220  FORMAT!////) 

1230  F0RMAT(I4,2X.F9.1,2X.F10.2) 

1240  FORMAT  ("SURGE.  WAVE  HEIGHT  OR  FALL  VELOCITY  BELOW./. 

$  'MINIMUM  ACCEPTABLE,  RE-ENTER  DATA*) 

1250  FORMAT! "TEMPERATURE  NOT  ACCEPTABLE .  RE-ENTER  DATA*) 

1260  FORMAT! "OUTPUT  RESULTS  TO  FILE  OR  LINE  PRINTER?  -  Y  OR  N*) 

1270  FORMAT!//. IX, "ESTIMATED  EROSION  ADJUSTED  FOR". 

$“  DURATION”. /.5X.*! IF  SUFFICIENT  MATERIAL  EXISTS)*,/, 

*/, IX. "EROSION  ABOVE  CHANGE  DATUM  (AS)  =  FR0M*,F0.2. 

<"  TO  \F3.2,1X.A2."3/*.A2,//) 

RETURN 

410  CLOSEiS) 

0PEN<5.FILE= ' INPUT ' ) 

RETURN 

END 

C 

C 

SUBROUTINE  SPEED 

REAL  PDIA , CTENP . FALL VEL . RHOSED, RHOSEA, K VCSEA, HU.RHD. K VC , GRNBUOY . R 
COMMON  /FALVZL/PDIA , CTEMP . FALLVEL . SURGE . HAVEHT.tEMP , BAREA 
DATA  RHOSED . RHOSEA , KVCSEA/2 . 70 . 1 . 027 . 1 . 048/ 

C 

C  GIVEN  MEDIAN  PARTICLE  DIAMETER  (IN  MM)  AND  HATER  TEMPERATURE 
C  (IN  DEG.  0  'SPEED'  COMPUTES  THE  PARTICLE  FALL  VELOCITY 
C  (IN  METERS/SECOND).  PROGRAM  ASSUMES  SEDIMENT  DENSITY 

C  (RHOSED  :  2.70  GM/CU.  CM).  AMBIENT  HATER  DENSITY  RATIO 

C  (SALTWATER.  33  PPT  -  RHOSEA  -  1.027)  AND  AMBIENT  HATER 

C  VISCOSITY  RATIO  (KVCSEA  =  1.048).  REFER  TO  CETN-II-4 

C  TOR  FURTHER  INFORMATION. 

ML1-. 0020319+. 015883/EXP  (CTEMP»*. 9/22. ) 

RHO--R!  1032 A»  ( 1  - .  0000 1 9547* ADS  (CTEMP  4 . ) ** 1 .68 ) 

KVC-KVCSEA*MU/RHO 

GRNBUOY = . 9806* ( RHOSED-RHO ) *PDI A**3/ ( RHO*KVC**Z ) 

R'K VC/PD IA 

IF ( GRNBUOY. LT. 39.) THEN 
FALLVEl=GRNBU0Y*R/130. 

ELSE  IF ■ GRNEUOY ,LE. 10000. )THEN 
rALLVELsGRNBU0Y**,7*R/60. 

ELSE 

FAL  LVEL-SQRT ( GRNBUOY ) tR/7.54 

END  IF 

RETURN 

END 

V. 

C 

SUBROUTINE  EXPLAIN 
C 

C  THIS  SUBROUTINE  EXPLAINS  THE  VARIABLES  USED  BY  THE  PREDICT 
C  PROCEDURE. 

C 

COHMON/BLK/IDATE. IT  ME. IPROF, ILOC. ISUR.DIST, IUNIT 
C 

C  DETERMINE  JF  HATER  TEMP  TO  BE  INPUT  IN  DEGREES  C.  OR  F. 

C 

ITEMP-1HF 


ir(IUNIT.EQ.2HN  HTEMPMHC 


C 

UR  JTEtt,.  100>  IUNIT.  IUMIT.  IUMIT 
IOC  FORMAT i 1X,/.5X.‘THE  PROCEDURE  REQUIRES  THE". 

$"  FOLLOWING  INPUT",/, 

$<•■/  " -  - -  -  —  -  —  - "  // 

12X’.  "SURGE  DATUM  -  THE  VERTICAL  DATUM  USED  TO  COflPUtE  THE  SURGE" 

*, /.8X. “LEVEL  (IN  " ,A2,"  RELATIVE  TO  THE  PROFILE  DATA)”,//, 

12X. "SURGE  HEIGHT  -  SURGE  HEIGHT  ABOVE". 

$"  THE  SURGE  DATUM  (IN  *.A2. “).”,//. 

$2X, "WAVE  HEIGHT  -  SIGNIFICANT  DEEPWATER  WAVE  HEIGHT  ", 

$' ( IN  *.A2. ").”./) 

C 

WRITE(t,110)IUNIT. ITEHP. IUNIT 
110  FORMAT i 

»2X. "CHANGE  DATUM  •  VERTICAL  DATUM.  (IN  ",A2. 

$'  RELATIVE  TO  THE  PROFILE  DATA)"./. 

$8X."  ABOVE  WHICH  THE  ERODED  AREA  WILL  BE  COMPUTED".//, 

$2K. "MEDIAN  SAND  SIZE  -  ALWAYS  IN  MM.*.//. 

I2X. "WATER  TEMPERATURE  -  IN  DEGREES  ",A1 ,//, 

$2X . "’ F ALL  VELOCITY  -  IN  *.A2."/SEC.  (OPTIONAL.  “ 

I, "PROGRAM  WILL  COMPUTE  THIS)”,/) 

WRI7E(&.115) 

115  FORMAT !2X, "STORM  DURATION  HOURS  SURGE  EXCEEDED  A  LEVEL"/. 

18X.“3.3  FT  (1  M)  LESS  THAN  THE  PEAK  SURGE."./! 

C 

WR!TE(6.117) 

117  FORMAT!  "CONTINUE?1') 

READ  ( 5 . 1 1 0 .  END  -1 35 .  ERR-- 1 35 )  IREPLY 

118  FORMAT ( 1A1 ) 

IF»  IREPLY.  NE.ltlY!GOT0140 
C 

HRlTEii.120) 

120  rORMAT!/.5X."THE  PROGRAM  OUTPUTS  THE  FOLLOWING"./, 

■ _ _ . _  _  ... _ _ _ ■  II 

12X,’ "TOTAL  EROSION  (AE)  THIS  IS  THE  VOLUM^  OF  THE*./ 

*.2i(.’  TOTAL  ERODED  AREA.  REGARDLESS  OF  DATUM".//, 

A2X. "ABOVE  DATUM  CHANGE  (AS)  -  THIS  IS  THE  VOLUME",/, 

J2X , "  CHANGE  ABOVE  THE  CHANGE  DATUM .",//. 

$ZX. "HORIZONTAL  SHIFT  -  THIS  IS  THE  HORIZONTAL"./, 

JZX."  SHIFT  OF  THE  CHANGE  DATUM", 
i"  INTERCEPT"./,/. 

12X.T0SITI0N  Or  THE  SURGE  INTERCEPT  ••  THE  POSITION."./, 

1ZX."  RELATIVE  TO  THE  SURVEY  DATA.  OF  THE"./, 
tZX,"  POST-STORM  SURGE  LEVEL  INTERCEPT. "./) 

C 

WEI7E(6, 130) 

130  FORMAT (2X. "ADJUSTED  EROSION  AMOUNTS  THIS  IS  THE  RANGE  OF",/, 

$8X. " -AS-  ADJUSTED  FOR  THE  DURATION  OF  THE  STORM  (AS  DEFINED1. /,8X, 
t” ABOVE).  MAXIMUM  ALLOWED  ADJUSTMENT  IS  50  PERCENT  OF  -AS-.".//) 
140  RETURN 
135  K7=E0T (5) 

RETURN 

END 
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APPENDIX  B:  NOTATION 


B1 


A  Volume  of  erosion  and  accretion  zones  computed  between  prestorm  and 
e  poststorm  profiles  (1/2  A^) 

A  Gross  change  between  prestorm  and  poststorm  profiles 
g 

A  ,  Volume  of  material  eroded  above  mean  sea  level  (msl) 
msl 

Ag  Volume  of  material  removed  above  the  surge  level 
Median  sediment  size 
Significant  deepwater  wave  height 
Lq  Significant  deepwater  wave  length 

Predicted  poststorm  dune  slope 
m^  Seaward  slope  of  the  poststorm  profile 
R  Maximum  shoreward  extent  of  poststorm  profile/surge  level  intercept 
S  Peak  surge  level  relative  to  mean  sea  level 
t  Water  temperature  during  storm  event 

w  Fall  velocity  of  median  sediment  size, 
x  Horizontal  position,  positive  in  seaward  direction 
y  Vertical  position,  positive  with  increasing  elevation 
a  Standard  deviation  of  prediction  of  A 

A  S 


B2 


END 

FILMED 
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